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Abstract 

r~| . Entanglement, which is an essential characteristic of quantum mechanics, is the key element in 

potential practical quantum information and quantum communication systems. However, there 
C^ ■ are many open and fundamental questions (relating to entanglement measures, sudden death, 

O^. etc.) that require a deeper understanding. Thus, we are motivated to investigate a simple but 

non-trivial correlated two-body continuous variable system in the absence of a heat bath, which 

> 

l/^ ' facilitates an exact measure of the entanglement at all times. In particular, we find that the 

^' 

On ' results obtained from all well-known existing entanglement measures agree with each other but 

^H ■ 

iy-\ , that, in practice, some are more straightforward to use than others. 

o 
o 



^ Dedicated to the memory of Krzysztof Wodkiewicz. 
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I. INTRODUCTION 

Entanglement, which is an essential characteristic of quantum mechanics, plays a key role 
in all applications related to information science [1-4]. Entanglement describes correlations 
between two or more particles or subsystems (qubits, oscillators, etc.). Despite the fact that 
much insight has already been obtained, much remains to be done as underlined by the title 
of a recent book on the subject [5]. For example, whereas there is a plethora of entanglement 
measures, there is also a consensus that a more unifying fundamental measure needs to be 
developed. Also, the discovery of "entanglement sudden death" (ESD) [6-8], in contrast to 
the well-known exponential decay of decoherence, requires a deeper understanding. Since 
investigations of ESD have also incorporated heat bath effects, master equations have been 
the tool of choice, despite their inherent limitations [9, 10]. Thus, we are motivated to analyze 
a simple but non-trivial correlated system which displays entanglement in the absence of 
a heat bath. For this system, an exact measure of entanglement exists which provides us 
with a touchstone for judging some of the various entanglement measures discussed in the 
literature. Since some of these measures involve entropy considerations, this should also 
throw some light on whether there is a close relationship between the entanglement of a 
system and its entropy. 

Thus, in Sec. II, we will consider a system of two free particles in an initially entangled 
state amd we calculate its time dependence. In Sections III, IV and V, we examine the 
same state using various entanglement criteria and demonstrate explicitely the various steps 
needed to demonstrate entanglement at all times. In Sec. VI, we present our conclusions 
that the various entanglement criteria lead to the same results but that the logarithmic 
negativity entanglement criterion is the simplest to use. 

II. ENTANGLED WAVE FUNCTION 

We consider two free particles, each of mass m, at positions Xi and X2, in an initially 
entangled Gaussian state. Thus, we are dealing with a system with continuous degrees 
of freedom (as distinct from a system of discrete variables such as qubits), applicable to 
particle position or momenta or to the field modes of light (of interest in connection with 
hnear optical quantum computing). 



The most general initial Gaussian wave function is 

, I ^^ (aiia22 - 'A-i)^^^ f «iia^i + 2012X1X2 + 022^2 1 .„ , s 

^(a;i,X2;0) = -= exp <^ ^ \. (2.1) 

In order that this state be square-integrable we must of course assume that an and 022 are 
positive and that 011022 ~ (^\2 > 0- We specialize to the symmetric case, by choosing 

an = 022 = ^ + ^, «-=(-^ + 4^)<0' (2-2) 

so that 2d > a (the opposite choice would lead to similar conclusions). As we shall see 
presently, d corresponds to the width of the center-of -mass system. Thus, 012 provides a 
measure of entanglement at t = 0. In fact, since we are dealing with free particles we expect 
that the entanglement will not change in time. However, as we shall see, the time-dependent 
coefficient of a;ia;2 will not serve as a measure of entanglement. 
The above equations enable us to write (2.1) in the form 

Next, we transform to center of mass and relative coordinates [11] 

a; = a;i-a;2; X = {xi+X2)/2 (2.4) 



m^ = m/2; M = 2m (2.5) 



to give the result 



^.fe...;O) = 0W*(.Y) = ^exp{-^-^}^ (2.6) 

It is clear that a is the width of the relative coordinate system whereas d is the width of the 
center-of-mass coordinate system. 

For future reference, we also note that the corresponding relative and center of mass 
momenta are 

p=7;iPi-P2); P = Pi+P2, (2.7) 



respectively. It is clear that both {x,p) and {X,P) satisfy the usual commutation relations 
for conjugate canonical variables. 

Since our wave function in the transformed coordinates now behaves as the product of 
two independent Gaussian wave packets, we can now apply the exact propagation method 
for a free particle [11] to obtain the wave function at a time t, with the result 
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where 



and 



It follows that [12] 



a\t) =aUl + 
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(2.9) 



(2.10) 



(x(t)) = 0; (X(t)) = 0; (M^)) = 0; (P(t)) = 



(2.11) 



{x'{t)) = a\ty, {X\t)) = al{t) 



(2.12) 



and 



We also note that 



<^^(^))-^^ <^^(^))-^- 



(2.13) 



aii(t) = a22(t) 



a2(t) AaM 



(2.14) 



In addition, we denote the coefficient of the real part of the (— |a;ia;2) power in the expo- 
nential in (2.8) by ai2{t) to obtain 



However, it is now clear that ai2{t) can not be used as a measure of entanglement because 
it predicts that entanglement decreases in time, eventually falling to zero [see the following 
section, especially (3.16)] whereas we expect that entanglement does not change in time. 
On the other hand, we have the tools to investigate the results obtained by use of various 
entanglement measures discussed in the literature. 

More generally, we point out that the starting point for all investigations is the well-known 
separability condition 

P = J2^d^^pf^ (2.16) 

j 
where p is the density matrix of the quantum state which is written as a convex combination 

of tensor product states for the j states and where the individual terms are normalized so 

that J2j Cj = 1- If this decomposition is not possible, then we say that the state is entangled. 

We now turn to some specific entanglement criteria which have been proposed. 

III. DUAN ET AL. [13] CRITERION 

Using the uncertainty principle, Duan et al. [13] derived a sufficient criterion for insepa- 
rability for a pair of EPR type operators for continuous variable systems. Further work on 
this topic appears in [14, 15]. In particular, this approach applies to our problem. Following 
[13], we write 

u = — l\a\xi + -X2J , V = — (\a\pi p2J. (3.1) 

(except that we have introduced the parameter L which has the dimension of length), where 
a is an arbitrary, non-zero real number. Use of the uncertainty relation [13] leads to the 
result 



±(A.V^(A"^>>(«^4 



(3.2) 



for any L. For L = 1, this is the Duan et ah result. However, by minimizing with respect 
to a and L, we will obtain an improved version given by (3.11) below. 
For a state that is not entangled, form 

{u') + {v^) = ^M) + ^(^?) + ^M) + ^M)- (3-3) 

Here we have used the fact that for a non-entangled state (2.16) the quantities 



(xixa) = (a;i)(a;2) = 0, (3.4) 

{Pm) = (Pi)(p2) = 0, (3.5) 

because we have restricted our discussion to states for which 

(xi) = (0:2) = (pi) = (P2) = 0. (3.6) 

(Duan et al. do not make this restriction but come to the same conclusion for what they 
call Am and Af .) Next we use the uncertainty principle, 

to get 

Now the quantity |/ + ^,0<t/<oo has a minimum value of 1 at y = |. We conclude 



{u") + (t;2) > a^ + — , (3.9) 



1 

independent of L. This is the result of Duan et al. Since the inequality is independent of L, 
we can minimize the left hand side with respect to L to get 



A 




(3.10) 



This is our improved inequality. It is a sufficient condition that the state is separable (not 
entangled). If it fails, the state must be entangled. We note that a^ + -^ has a minimum 
value of 2 at a = ±1. Hence 
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^ h. (3.11) 



This is the necessary condition that a two-particle state be separable. Thus, we have two 
necessary conditions, corresponding to choosing a to be positive or negative. Using (2.4) 
and (2.7), we may write the two conditions in the succinct forms 

(^^)(P^)^^(^), (3.12) 

and 

(x^)(P^)^4(^). (3.13) 

We now with to apply these results to the particular state discussed in Sec. II. Thus, using 
the results given in (2.13) and (2.15), together with (2.9) and (2.10), these conditions take 
the explicit forms 



and 



.^+f^)%«^. (3,15) 

\m(j J 

Since we assumed 2d > a, it follows that (3.14) (which corresponds to the choice of positive 
a) is automatically fulfilled, implying separability. However, (3.15) (which corresponds to 
the choice of negative a) is only fulfilled if 



f2mad\ f a^V'" , , 



1/2 
1/2 



Thus, for t < tfii the separability condition is violated and the state is entangled. However 
at t = trf, we encounter ESD [6-8], despite the fact that we know from our exact analysis in 
Sec. HI that the state is entangled for all times. The solution to this apparent contradiction 
stems from the fact that the Duan et al. condition is a sufficient criterion for inseparability 



(entanglement) but it is not necessary. Recognizing this, these authors were led to develop a 
necessary and sufficient condition for entanglement by using a variety of local linear unitary 
transformations (consisting of various rotations and squeezing transformations of Xi and pi 
that preserve the commutation relations and similarly for X2 and P2), and referred to as 
LOCC, to map any Gaussian state into what they refer to as Standard forms I and II, which 
eventually leads them to a state for which their separability criterion is both necessary and 
sufficient. It is known that these local operations do not affect the entanglement of the state 
i.e. we have a family of states all with the same entanglement [4, 16]. 

Guided by the fact that our results are satisfactory at t = 0, supplemental by the detailed 
results which we already obtained for the motion of a free particle [12] , especially equation 
(9) of the latter reference, we make the following local canonical transformations: 



Xi{t) = exp{--Hit)xi{t)exp{-Hit); X2{t) = exp{--H2t)x2{t) exp{-H2t), (3.17) 



where 



Hi = pl/2m; H2 = pi/ 2m. (3.18) 

It follows that [since pi(t) = pi{0) and P2(t) = p2( 



xi{t)=xi{t)--pi{0)=xi{0), (3.19) 

m 

X2{t)=X2it)--p2{0)=X2{0). (3.20) 

m 
In fact, the second equality in the latter two equations readily follows from the Heisenberg 

equation of motion. These transformations lead to the results 

x(t) = x(t) p(0) = x(0), (3.21) 

X(t)=X(t)-^P(0) = X(0). (3.22) 

It follows that the necessary conditions for separability now becomes 

{X'){p') > \ (^) (3.23) 



and 



{x'){P')>A{^] (3.24) 



which takes the exphcit forms 



2d>(T (3.25) 



and 



(J > 2d, (3.26) 

which are only compatible for a = 2d. However, since we assumed that 2d > a, it is clear 
that the above analysis leads to the conclusion that the system is entangled for all times. 
We note that the unitary transformations given in (3.17) led to the elimination of terms 
depending on t in the separability conditions. 

A similar analysis may be carried out using the Peres-Horodecki criterion [17, 18], namely 
that a state is separable if the partial transpose of the density matrix is a positive operator. 
It can be shown explicitly (See Appendix A) that, for our model, it leads to the same result 
(3.11), that was obtained from the Duan et al. criterion. 

Although the model we are considering here is a pure state, it is instructive to see how 
it fits into the general framework of mixed states which are best considered using Wigner 
distributions [13, 21]. 

It should be emphasized at the outset that not all Wigner functions are permissible 
distribution functions since the corresponding density matrix elements must be positive 
definite [22] and the uncertainty relations must be satisfied. We now briefly review the work 
of Duan et al. [13] and Simon [21], which will result in bringing M into the "standard form" 
[13] given in (3.30) below. 

Recalling that Gaussian states are completely characterized by their first and second 
moments (and here we have arranged that the former are zero), it follows that the Wigner 
characteristic function for a Gaussian state of a pair of particles can be written in the general 
form 

W{Qi,PuQ2,P2;t) = exp{- ^'^'^ ], (3.27) 



where 



ft 



Q 



ft 



M 



G C 

C^ H 

Here M is the correlation (variance) matrix and G and C are 2x2 matrices given by 



(3.28) 
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H 



2ft 
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2ft 



2ft 

ft2 

(3;2P2+P23:2) 
2ft 



'W 



L2 ft 

(3:^2^1) L'^{pip2) 
ft ft2 



(3.29) 



In these expressions L and h are constants introduced to make the matrix variance 
(correlation) M dimensionless. However, as far as the subsequent analysis is concerned, 
the L may be ignored with impunity, as we will do henceforth. 

Making use of a series of local linear canonical transformations (rotations and squeezings), 
it was shown [13, 21] that it is possible to bring M to the special form: 



M' 



(3.30) 



f g c 0\ 

g c' 

c h 

\0 d h J 

Since determinants are invariant under these transformations we have the following simple 
relations for determining the quantities g, h, c and c', in terms of four invariants. 



detG 

detC 



9 , 
cd . 



det H = /i^, 
det yi= {gh- c^) [gh - c'^) . 



(3.31) 



We now turn to the special case of interest here, that is the pure Gaussian state given in 
(2.1). Using the techniques developed in [23], or, since we are dealing with a free particle. 
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from the Wigner function given in (A5), together with the results given in (2.14) and (2.15), 
it follows that 



Hence 
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It follows that 
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(3.32) 



(3.33) 



(3.34) 



det G = det H = ( — + — r 1 a^d 
a'^ Ad^ 



2j2 



4 V 0-2 4^2 



(3.35) 
(3.36) 



det M = det(G + C) det(G - C). 



(3.37) 



We note that the latter two equations are independent of t. Thus, the transformed matrix 
M has the form (3.30) with 



1/1 1 



^" = 4^ + 4^' "^ 



2j2 



(3.38) 
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The solution of these equations is 



-'-iB+4^)^^^'^ ^'-''^ 



(/ - c^) (/ - C-) = 4- (3-40) 



^=r^(+^+i) ^'-''^ 



c = -c= -lad (-^ + 7^ ) • (3.42) 



In terms of these quantities, the inequahty (3.11) becomes 



^{g Tc){g± d) ^ 1. (3.43) 

In terms of the above expressions, this becomes 

which imphes 

a ^ 2d, or 2d ^ a. (3.45) 

In other words, the condition for separability only holds when 

2d = a (3.46) 

that is when the center-of-mass coordinate width and the relative coordinate width are equal. 
As a consequence, 

ai2 = 0, (3.47) 

as expected. As indicated above, these conclusions hold at all times. 

It is also of interest to note that, within the present context, the Peres- Horodecki criterion 
implies that a Gaussian state is separable if and only if the minimum value of its symplectic 
spectrum of M^^'^ is greater than 1/2 [19, 21] which leads to a good measure of entanglement 
for all Gaussian states given by 

12 



E = max {0, - log(2z/™„)} (3.48) 

where Umin is the smallest sympletic eigenvalue of M'^'^ . The equation determining the 
sympletic spectrum is [19] 

z/^ + (det G + det H - 2 det C) u^ + det M = 0, (3.49) 

with solutions ±iz/Q,, a = 1,2 where u^ is the symplectric spectrum. Hence, using (3.35) 
and (3.36), we obtain z/i = (d/a) and 1^2 = {a/4d). Since we assume 2d > a, we see that 
^min = ^2 < \ a^d heuce 

E = max|o,log( — J I =log( — J , (3.50) 

in agreement with the result (4.10) arising from the log negativity criterion, as discussed in 
the next section. 

IV. LOGARITHMIC NEGATIVITY CRITERION 

The logarithmic negativity is defined as 

i?^(p) = log{2iV(p) + l}, (4.1) 

where N{p) is the negativity of the state and is given by the absolute sum of the negative 
eigenvalues of the partial transpose of p [19, 20]. 
We want to solve the eigenfunction equation: 

dxi dX2{x[,X2\p'^^\Xi,X2)(j){Xi,X2) = \(l){x[,X2), (4.2) 

especially in order to obtain the negative eigenvalues. After some algebra (See Appendix 
B), we find that the eigenvalues are given by 



±/3"+" Ao for m ^ n 
\mn = { (4.3) 

/3^" Ao for m = n, 

where m, n = 0, 1, 2, and Aq is a positive eigenvalue given by 
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In addition, 



/3 = W'^^^^^7==^- (4-5) 

flu + V Oil ~ '^12 



For the symmetric case, these reduce to 



" (2rf + a)2' 
and 



recalhng that we have assumed that 2d > a. As a check, we note that 

Y, Amn ={l + (5^ + (5' + ..)\,= (y^) Ao = 1, 
verifying that Trp^"^ = 1. In addition, 



Hence 



(4.6) 



^-1^- (-) 



(4i 



(4.9) 



En{p) = log{2iV(p) + l} 

= log 1^1. (4.10) 

Thus, the greater 2d is compared to o", the larger the negativity and hence the greater the en- 
tanglement. In addition, since p{t) = exp {—iHt) p(0) exp (iHt), where H = {pi + P2) /2m, 
it is clear, from (4.2), that the eigenvalues of p^^ are invariant under this local unitary 
transformation. Hence, the result (4.10) is valid for all times. 
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V. ENTANGLEMENT OF FORMATION 



For bipartite pure states, the entanglement of formation is given by [24] 



Ep = 5*1 (pi) = 5*2 (P2) 
where pi and p2 are the reduced density matrices [defined in (B15)] and 



(5.1) 



Si = -Tr2[pilogpi] 
is the von Neuman entropy. Thus, from appendix B, we have 



(5.2) 



S = -^XJog \r, 



n=Q 



where the eigenvalues of pi are given by 



A„ = Ao/3'" (n = 0,l,2, ), 

where, from (4.6) and (4.7), we have (with R = 2d /a) 



and 



Hence 
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(5.5) 



(5.6) 



^1 = -AoX]/3'" [logAo + 2nlog/3] 



n=0 



- Ao log Ao Y, /?'" - 2Ao log [5 Y, <P 



2\n 



n=0 



n=0 



2/32 

- log ^0 - 7^ log /3 



logi? + 



1-/32 

'^ l + R iR-l)\ R + l 
2 log ^77^ + ^„ log 



(5.7) 



2R 2R ° R-l 

We recall, from (4.10), that the logarithmic negativity En{p) is given by log-R. Also it 
can be shown that < 5*1 < logi? for i? > 1 and the equalities hold for i? = 1 and i? — )■ 00, 
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respectively. This is consistent with the result that the entanglement of formation is always 
less than logarithmic negativity, and they are equal for maximally entangled pure states. 
Thus, the entanglement exists if 2d > a. 

VI. CONCLUSIONS 

We examined a simple but non-trivial model of entanglement which enabled us to carry 
out an exact analysis. We analyzed various entanglement criteria, arising especially from the 
work of Duan et al. [13]; Peres-Horodecki [17, 18]; Vidal and Werner [19], who considered 
both the logarithmic negativity and that arising from a determination of the smallest sym- 
pletic eigenvalue of the Peres transform of the transformed variance matrix and Bennett et 
al. [25] on the entanglement of formation. We found that all of these various entanglement 
criteria led to the same results but that some are more straightforward than others. In par- 
ticular, it was clear that the logarithmic criterion is the simplest to use since the procedure is 
straightforward, that is obtain the eigenvalues of the Peres transform of the density matrix. 

After this paper was completed, we became aware (courtesy of the referee) of various 
papers that have closely related themes. Our work is an example of "entanglement without 
dissipation," which apparently was initially discussed by Chan and Eberly [25] who also 
investigated a Gaussian state but used a Schmidt-state analysis as a measure of entangle- 
ment. Next, Yonac et al. [26] considered two isolated atoms each in their own lossless 
Jaynes-Cummings cavity. They showed that, due to the interaction with the local lossless 
cavities, ESD occurs for atom-atom entanglement due to information loss to the cavity modes 
but that entanglement is resurrected in a periodic manner following each ESD event due to 
the fact that the time evolution is lossless. The same system was analyzed by Sainz and 
Bjork [27] who concluded that the atoms simply transfer their entanglement to the cavity 
fields and that an entanglement measure exists that is constant under the time evolution. A 
different system, photoionization in a lossless environment, was considered by Fedorov et al. 
[28], who found narrowing of electron and ion wave packets due to electron-ion entanglement. 

Entanglement of formation [24] is one of the measures we have discussed (see Sec. V) 
and this quantity is referred to by Munro et al. [29] as "- - the canonical measure of 
entanglement - -," who then go on to present a class of states that have the maximum 
amount of entanglement for a given linear entropy. All of this work is leading to a better 
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understanding but, to quote from the recent general overview of Yu and Eberly [30], "- - 
there is still no deep understanding of sudden death dynamics." 
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Appendix A: Density Matrix Elements and Wigner Functions 

In general, the Peres partial transpose of the density matrix is 

(x'-^^x'2\p^'^\xi,X2) = {x[,X2\p\Xi,X2) ■ (Al) 

However, in practice, it is often more convenient to consider the corresponding result for the 
Wigner function [21], that is 

^^^'{ql,Pl■,q2,P2) = W{ql,Pl■,q2,-p2). (A2) 

The corresponding transpose of the Wigner Characteristic function is 

W^'{Qi,Pi; Q2, P2) = W{Qi, Pi; -Qa, P2), (A3) 

We recall that the most general Gaussian pure state corresponds to the wave function 

(011022 - 012)^^"^ ^^^ / aiixj + 2ax2X\X2 + a22x\ \ 



i,{x,,X2) = '^''^'' ^''' cxp <i - --^^ — --- — ^ 1 , (A4) 



where a\\ and 022 are positive and an 022 — O12 > (to ensure integrability) . It is straight- 
forward to obtain the Wigner function 



..,. N 1 / oiig^ + 2ai2gig2 + 022^1 1 
>V(gi,Pi;g2,P2) = T^j^^^P] ^ f 

xexp|-2^H4^1^^iHEl^i±^\ (A5) 

1^ n^ (011022 - ofs) J 

and the Wigner characteristic function (Fourier transform of the Wigner function) 
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W{Qi,Pi;Q2,P2) = exp 



1 O22-P1 — 2ai2PlP2 + ^11-^2 I / \ n\ 

X exp <^ '- ^-— '- } . A6 

2{aiia22 - a{^)h? 



Also, the corresponding density matrix is 



(x'l, X2|p|Xi, X2) = i) K, X2) ^*(Xi, X2) 
_ (^11022 — ^12) 

/ an (x^ + xf ) + 2ai2 (xia:2 + x[x'^) + 022 (a:| + xf ) . 
The Peres partial transpose of this density matrix is 



/ / / I T2 1 \ III \ I* I l\ ('^ll'^22 — ^12) 
{X^,X2\P \Xi, X2) = i! (Xi, X2) i! (Xi, X2) = 

f au(xl + xf) + 2ai2(xix'2 + x[x2) + a22(xl + Xn)'] , . ^. 

xexpj |. (A8) 

The corresponding transpose of the Wigner characteristic function is 



W^HQ,,P,;Q2,P2) = exp|-^H^i^^^iH|i^i±^HH^| 

X exp{-^H4f^^iH^^4±^}. (A9) 

In addition, the corresponding transpose of the Wigner function is 



>V^'(gi, Pi; 92,^2) = — -2exp<^ 
{-Kh) {, 



\_ „„„ [ '^11^1 + 2ai2gig2 + 022^1 



xexp|-2 "^^l+'""^^^^t!"^H - (AlO) 

l^ n^ (011022 - afs) J 



Consider the symmetric case, for which 



«- = «n = ^ + i^' «- = -^ + i^- (All) 



Then 
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{x[,X2\p ^|a;i,a;2) 



27r 



'2\ 



, aii(a:^ + xf) + 2ai2(a:ia:2 + 3:13:2) + 011(3:2 + x'2 j , ... 



iy^^(gi,Pi;Q2,P2) = exp 



iQi + Q2f (Q1-Q2)' 



5(7 



2 



32rf2 



72 ^2 



X exp{-^ (Pi + P2)' - ^ (^1 - P^n (A13) 



>V '(9i,Pi;?2,P2) = -—2 exp 



{nhf ^ \ 2^2 8rf2 

X exp <( --- (pi - p^f - 7^ (pi + P2)' I . (A14) 



^2 v^^ ^-/ 2^12 

If W^^ is to be a Wigner function we must require that the uncertainty relation be 
satisfied. In particular we must require 

AxAp > -, (A15) 

where 

{Axf = {{xi -Xif) ' 

dqi / dpi / dq2 / c?P2 (^i - ^2)^ >V^'(gi,Pi; ^2,^2) 

dqi / rfpi / rfg2 / dp2{qi - ^2)^ W (gi,pi; ^2,^2) 
. = ((3:1 -3:2)^), 
A{Apf = {ipi-p2?f' 

dqi / dpi / rfg2 / c?P2 (Pi - ^2)^ >V^'(gi,Pi; ^2,^2) 

dgi / dpi / (ig2 / c?P2 (pi + ^2)^ W (gi,pi; ^2,^2) 
= ((Pi+P2)'>. (A16) 

Therefore, if we require that the Peres transpose corresponds to a Wigner function, we see 
from (A15) and (A16) that the following inequality must be satisfied 
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\J {{x, - x^f) {{p, + p^f) > h. (A17) 

This is just the refined condition of Duan et al. 

Appendix B: Eigenvalues and Eigenfunctions of the Density Matrix and the Re- 
duced Density Matrix 

For the case of the pure Gaussian state given by (2.1), and speciahzing to the symmetric 
case, where 022 = an, we found that the Peres transpose is [see (A8)] 

{x[, X2 \p^^ I a;i, ^2) = ■?/' {x[, X2) ip* (xi, X2) , (Bl) 

where [see (2.1)] 

^(xi,4) = ("^^ -^^^'^ exp |-"" ^"^ + ""] + '"""^"^ } , (B2) 



with a corresponding result for ■il){x'i, X2) ■ We want to solve the eigenfunction equation: 

dxi I dx2{x'^,x'2\p^'^\xi,X2)^{xi,X2) = \^{x\,x'2). (B3) 

The form of (B2) suggests use of the Mehler formula [11], which is written in terms of the 
Hermite functions i/„. However, we find it is more useful to modify this formula so that it 
is now written in terms of the related orthogonal function 0„, the eigenstate of the quantum 
harmonic oscillator. Thus, for arbitrary variables x and y, the modified Mehler formula may 
be written in the form 



V L P L 1 ) „=o 



>n(a:)0„(y), (B4) 



where the function 0„(x) is related to the Hermite polynomial Hn{x) by 



and 



U-) = ^^ {D "'' e--V2^„ (^ ^) , (B5) 



(t)m{x)(l)n{x)dx = Smn- (B6) 
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Thus, applying (B4) to both ip^xi.x^) and ip{x'i,X2), we obtain 

(x;,x^|p^^|xi,a;2) = Ao5^/3"+>™(a;i)0^(4)0„(a;;)0„(a;2). (B7) 

Here the parameters 7, /3, and 70 are given by 



7 = sj< - alJ2 (B8) 



/^=J """^f "t (B9) 



Ao = ^^"7^ "'% =1-/3^. (BIO) 

On + V '^11 ^ '^12 

The advantage of writing the matrix elements in terms of orthogonal functions, as given in 
(B7), immediately leads to solutions of the eigenfunction equation (B3), in the form 



dxi I dx2{x[,x'2\p^^\xi,X2) 0m(a:i)0„(a;2) = Ao/3™+"0r«(a;i)0n(a;2)- (BH) 



If m = n, (^nn{xiiX2) = 0n(a;i)0„(a;2) is the desired eigenfunction with eigenvalue A„„ = 
/S^^Aq. On the other hand if m 7^ n, we find that the eigenfunction is given by 

$mn (a^l, 0:2) = -J= [(pm{Xi)(f)n{x2) ± (pn{xi)(f)m{x2)] (B12) 

v2 
with eigenvalues ±Ao/3™''^". In summary, the complete eigenvalues are 



i/J^+^Ao for m^n 

^mn = { (B13) 

/3^"Ao for 7n = n, 

with m,n = 0,1, . . . 

Next, we consider reduced density matrices. Thus, we define 

Pi=Tr2{x[,X2\p^''\xi,X2). (B14) 

This is obtained by setting X2 = x'^ in (B7) and integrating over X2 to get 

pi = Ao5^/3=^'"0™(xi)0„(x;). (B15) 
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Thus, the eigenvalue equation for the reduced density matrix is 

I dxipM^\xi) = A^^K). (B16) 

It immediately follows that this leads to solutions (f)m{xi) with corresponding eigenvalues 
Ao/3^™(^Ti = 0, 1,2 — — — ). Moreover, it is clear that the corresponding eigenvalues of p2 
are the same. 
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